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Abstract 

In this work, Emergent Universe scenario has been developed in general 
homogeneous anisotropic model and for the inhomogeneus LTB model. In the 
first case, it is assumed that the matter in the universe has two components 
- one is perfect fluid with barotropic equation of state p = ujp (to, a constant) 
and the other component is a real or phantom (or tachyonic) scalar field. In 
the second case, the universe is only filled with a perfect fluid and possibilities 
for the existence of emergent scenario has been examined. 
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1 Introduction 

The idea of emergent universe is the result of the search for singularity-free inflationary 
model in classical general relativity. An emergent universe model can be defined as a 
singularity free universe which is ever existing with an almost static nature in the infinite 
past (t — > — oo) and then evolves into an inflationary stage. In fact, an extension of 
the original Lemaitre-Eddington model can be termed as emergent universe. There are 
several features for the emergent universe viz. (i) the universe is almost static at the 
infinite past, (ii) there is no timelike singularity, (Hi) the universe is always large enough 
so that classical description of space time is adequate (iv) the universe has accelerated 
expansion, etc. 

However, there was singularity-free solution in the literature since 1967 - Harrison 
described a closed model of the universe with radiation, which coincides with Einstein 
static model in infinite past. But actual search for Emergent model of the universe was 
started ( about 40 years back ) by Ellis and collaborators (2004, 2004). They formulated 
a closed model of the universe filled with two non-interacting fluids - one is a minimally 
coupled scalar field having self-interaction potential and the other is a perfect fluid with 
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equation of state p = up. In fact,they studied only the asymptotic behaviour to charac- 
terize the emergent scenario without finding exact analytic solutions. Then Mukherjee 
ct. al. (2005) solved semi-classical equations in the Starobinsky model for fiat FRW 
space-time and examined the features of emergent scenario. Subsequently, Mukherjee 
and others (2006) were able to obtain nonsingular (i.e. geodesically complete) inflation- 
ary solution where a part of the matter is in exotic form. Mulryne et. al. (2005) have 
discussed the existence and stability of emergent models in the context of Loop Quantum 
Cosmology. Debnath (2008) has formulated an emergent model of the universe for exotic 
matter in the form of phantom or tachyonic field. Banerjee et. al. (2007, 2008) have 
shown a model of emergent universe in brane scenario while Campo et.al. (2007)have 
studied a model of emergent universe for self- interacting Brans-Dicke theory. Recently, 
Mukerji et.al. (2010, 2010) have formulated an emergent universe model in the context 
of Einstein-Gauss-Bonnet theory and in Horava Gravity. 

In section 2, we consider a homogeneous and anisotropic space time models described 
by the line element 

ds 2 = —dt 2 + a'dx 2 + b 2 dtl\ (1) 
where the scale factors a and b are functions of time t alone. We note that 

dfi 2 = dy 2 + dz 2 , when k = (Bianchi I model) 

= dd 2 + sin 2 6d(j) 2 , when k = + 1 (Kantowaski — Sachs model) 



= d6 2 + sinh 2 6d<f> 2 , when k = — 1 (Bianchi III model) 

Here k is the scalar curvature and the above three types are characterized by Thorne 
(1967) as flat, closed and open respectively. 

Now, the expression for the Hubble parameter H, the deceleration parameter q and 
the anisotropy scalar a in terms of the scale factors are : 



H 



a b 



q = -1 - 



and 



H 
IP 
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In section 3, we consider our universe as inhomogeneous Lemaitre — Tolman — 
Bondi[LTB] Model. The LTB model is one of the most wellknown sphercally symmetric 
models in general relativity. It was described by Lemaitre (1933), Tolman (1934) and 
Bondi (1947) during the period of time from 1933 to 1947. The exact solution have been 
obtained by Bonnor (1972, 1974). This simple inhomogeneous cosmological model agrees 
with current supernova and some other data (Moffat 2005; Moffat 2006; Alnes 2006; 
Alnes 2007; Romano 2010, 2010). Also very recently Clarkson and Marteens (2010) give 
a justification for inhomogeneous model from the point of view of perturbation analysis. 

The inhomogeneous spherically symmetric LTB space-time model is described by the 
metric ansatz in a co-moving frame as 

R' 2 

ds 2 = -dt 2 + — — — dr 2 + R 2 (d9 2 + sin 2 8d(f> 2 ) (2) 
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where R = R(r,t) is the area radius of the spherical surfaces and f(r)(> —1) is the 
curvature scalar that classifies the space-time as 

(i) bounded, if, -1 < f(r) < 0, 

(ii) marginally bounded, if, f(r) = 0, 

(iii) unbounded if f(r) > 0. 

The energy-momentum tensor for perfect fluid is given by 

T M „ = (p + p) u^u v + pg^ (3) 

where the fluid 4- velocity u M has normalization u M u M = — 1 and p, p are respectively 
the usual matter density and pressure of the fluid. Finally, the paper ends with a brief 
discussion in section 4. 



2 General anisotropic model : Basic Equations 
for Emergent Scenario : 

The Einstein field equations for homogeneous and anisotropic model of the universe with 
non-interacting two fluid system can be written as : 

^ + 2^ = — ^ [p m + p,f, + 3p m + 3p^] (4) 

and 

b 2 n ab k , , 

¥ + 2 ab + ¥ = [p - + P * ] (5) 
where p m and p m are the energy density and pressure of a perfect fluid having equation 
of state pm = Lopm, (w,a constant), p^ and pj, are the energy-density and pressure of a 
scalar field <f> having expressions 

P4, = 1 4> 2 + V(4>) (6) 
= |<^ 2 - Vtf) (7) 

Here V(4>) is the potential for the scalar field 4> and e — ±1 corresponds to normal or 
phantom scalar field. 

As the two components of the matter field are non-interacting, so the energy conser- 
vation equations are given by 

p m + 3H{p m +p m ) = (8) 

and 

p^ + 3H(p^+p^) = (9) 

For emergent scenario the appropriate form of the scale factors will be according to 
Debnath (2008) 

o(t) = oo [0 + e at ] n (10) 

and 

b(t) = b [5 + e !t ] m (11) 

where do , /3 ,a , n and &o , 8 ,7 , m are positive constants. The justification for such 
choices are the following : 
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1. ao > and bo > for scale factor to be positive definite. 

2. p > and S > 0, otherwise there will be big rip singularity. 

3. a > 0, n > 0, 7 > and m > for expanding model of the universe. 

4. a < 0, n < or 7 < 0, m < implies that there was a big bang singularity at 
infinite past. 

For the above choice of the scale factors, the universe started with a finite volume at 
t = — 00, grows gradually without encountering any singularity for any t and finally, the 
universe will be of infinite volume at future infinity. 

For this choice of the scale factors 'a' and 'b' , the Hubble parameter, it's derivatives 
and the deceleration parameter ' q' have the following expressions : 
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H = 



\P + e< 

{n(3(a 2 )e° 
[P + e at ] 



+ 



[2m 7 e 7t ] 
[8 + eT*] 

[2m5(7 2 ) 1 



3 



[nPa 3 e at (P - 



[P + e at f 



')] + [2m57~e 



[5 + eT*] 2 
3 ° 7 *(<5 



= -1 - 



3nPa 2 e at (5 + e 7 *) 2 + 6ro<Ve 7 *Q3 + e Qt ) 
[nae a '(5 + eT') + 2m7eT i (/3 + e Q ')] 2 



(12) 



We note that H and H are positive definite while q is negative definite throughout 
the evolution but H = when t satisfies 



(J + e 7t ) 3 Q3 - e at ) e 7 * _ 2m<? 7 3 
(jS + e at ) s (eT* - 5) e T* ~ n/3a 3 



Asymptotically, as t — > —00, ff, ii", and H all tend to zero but q — s> —00 while the model 
becomes a de Sitter universe as t —¥ 00 

In figure 1 we have plotted the deceleration parameter, q, against time ,t. As H is 
positive throughout the evolution so q is always less than —1 while asymptotically as 
t — > +00, H — > and hence from (|12l) . q — > —I as shown in figure 1. 

Now using the equation of state for the perfect fluid the energy conservation relation 
((8| can be integrated to obtain 



pm = po(ab ) 



-(l+w) 



(13) 



where po is an integration constant. 

Also, writing the explicit form of p$ and , from the field equations and © 



and 



i.e., we have 



P4> = 



IP 
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a b 



k_ 
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2a Ab b 2 2ab k 
a b b 2 ab b 2 



(14) 



(15) 



k_ 

V 2 



3H - 3o-" 



(1 + oj)p (ab 



(16) 
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Variation of Deceleration parameter against Time 
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t 

Fig.l 



Fig. 1 represents the graph of the deceleration parameter, q against the cosmic 
time,i, in anisotropic model of universe for the constants a — 0.1, /3 = 0.2, 7 = 
0.3, S = 0.4, m = 4, n = 6. 
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i(l-^po(a& 2 r (1+ ^ (17) 
it has the same value for real or phantom 



2.1 Possibility of Emergent Scenario: The restrictions. 

In this section, we discuss the possibility of emergent universe and present the restrictions 
for its validity. In other words, we examine whether the choices of a(t) and bit) given by 
equations (|10[) and (|11[) are possible solution for the field equations for (i)the real scalar 
field and (ii) the phantom scalar field. 



and 



v(4>) 



™ 2 + § 



We note that V(<j>) is independent of e i.e. 
scalar field. 



2.1.1 Real Scalar field (s = 

In this case, the expression for <j> 2 



-1) 



(1 + uj)p (ab ) 



-(!+") 



When the perfect fluid matter component is not of exotic nature (i.e. does not violate 
the weak energy condition) then for real scalar field, emergent scenario is not possible 
for k — 0, — l(i.e. for flat and open model) while for K ant awski- Sachs model (with 
k = +1), emergent scenaeio is possible provided 



o 2 . 3. . . 2\-C 1 + w ) 

3° + 2' + u)pa{ab ) 



(18) 



On the other hand, if the perfect fluid matter component is of phantom nature, then 
emergent scenario is possible for all three models provided 



\1 + oj\p Q {ab 2 ) (1+ " J) > 3H + 3a 2 



k_ 

b 2 



(19) 



and S can be obtained as 




b 2 2ab a 2b k 
b 2 ab a b b 2 



- {l + u)p {ab 2 )-^+")dt 



As a and b are given functions of t, (equation (llO|l and (|1 1 |l V so in principle the above 
integral can be evaluated to obtain as a function of t. 

The graphical representation of V for various choices of the parameters and the cur- 
vature scalar k are shown in figure fig 2(a)-2(h). The first five figures are for k — +1 then 
next two are for k = — 1 and the last one is for k — 0. The curves for the potential show 
some distinct features for different choices of the parameters. 



2.1.2 Phantom Scalar field (e 

In this case, d> 2 can be written as 



-1 



2 

3 L 



3H + 3a - — 
b 2 



(l + a;) /90 (ab 2 )- (1+ " ) 
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Fig.2a Fig2b 
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Fig.2e Fig2f Fig2g Fig.2h 

Fig. 2a represents the graph of the scalar field potential V against cosmic time, i, for 

the variables a = 0.05, (3 = 0.05, 7 = 0.1, S = 0.3, m = 1.05, n = 1.71, a = 
1, bo = 2, po = lj lj = 1/3 at k=l,i.e., for the closed universe. Fig. 2b represents 
the graph for the same value of the variables except for m = 1.38, n = 1. Fig. 
2c represents the graph for the same value of the variables except for m = 2.98. 
Fig. 2d represents the graph for the same value of the variables except for m = 
4.72, n = 1.71 Fig. 2e represents the graph for the same value of the variables 
except for m = 1, n= 2.48, a = 0.218, /3 = 0.228, 7 = 0.158, <5 = 0.3419,and 
a very low value of po = 0.001 Fig. 2f represents the graph of the scalar field 
potential V against cosmic time,i,for the variabless a = 0.812, (3 — 0.456, 7 = 
0.496, S = 0.596, m = 2.85, n = 2.95, at k=0,i.c., for the flat universe.Fig. 
2g represents the graph of the scalar field potential V against cosmic time, t, for 
the variables a = 0.608, /3 = 0.336, 7 = 0.1, 6 = 1.39, m = 2.81, n = 3.55, 
at k=-l,i-e., for the open universe.. Fig. 2h represents the graph of the scalar 
field potential V against cosmic time, i, for the same value of the variables except 
7 = 0.442 at k=-l,i.e., for the open universe. 
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when the perfect fluid matter component is not of exotic nature, the emergent scenario 
is possible for phantom scalar field when k = 0, — 1 (i.e. Bianchi I and Bianchi III model) 
while for Kantowski Sachs model (k=l) emergent scenario is possible provided 

1 < 3H + 3a 2 + |(1 + oj) Pi) (ab 2 y (1+ ^ (20) 

But when the perfect fluid matter component is of phantom nature, then emergent 
scenario is always possible for all three models provided 

\1 + aj\po(ab 2 y {1+hJ> <3H + 3a 2 - A (21) 

2.2 Emergent Scenario for Tachyonic Scalar field : 

For a tachyonic scalar field tp, with potential B(ip), the energy density p^, and pressure 
p^ are given by 

Py, - , (22) 
VI - eip 2 

and 

py, = -B(ipWl-siP 2 (23) 

where as before e = ±1 corresponds to normal and phantom tachyonic field. So, ip an d 
B(tp) can be written as 

i? = P -±±^± (24) 
e(PV-) 

and 

BWO = y/-Pi>Pi> (25) 
Hence from the filed equations ?/> 2 and B(ip) have the expressions 

= ,n j^-aH-^-d + o,)^)-^ 
H + ^ + ^-po(a6 2 )-^ 

and 



BW 



5b 2 ab o - fc 
2^ + ^ +3J/+ ^ 



+ o;po(afe 2 )-(i+") x ^ + — + ¥ - p (ab 2 y 

(27) 

We note that as before B(tp) is independent of e i.e. it has the same value for normal or 
phantom tachyonic field. 

The variation of B with cosmic time t for different choices of the parameters and the 
curvature scalar k are shown in fig. 3(a)-(h). The first four figures are for k — l,then two 
for fc = — 1 and the last two for k — 0. In this case also the curves for the potential show 
some special feature for different choices of the parameters. 



-(l+w) 
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Fig. 3a 



Fig3b 



Fig3c 



Fig.3d 




Fig.3e Fig3f Fig3g Fig.3h 

Fig. 3a represents the graph of the tachyonic field potential B against cosmic 
time, t, for the constants a = 0.1, (3 = 0.2, 7 = 0.1, 6 = 0.2, m = 1, n — 1, a,Q — 
0.2, bo — 0.2, po — 0.01, 10 = — 1 at k=l,i.e., for the closed universe. Fig. 3b 
represents the graph of the tachyonic field potential B against cosmic time, t, for 
the constants a = 1, (3 — 16.2, 7 = 0.59 6 — 0.2, m = 1, n = 1, ao = 0.2, fro = 
0.2, po = 0.01, = — 1 at k=l,i.e., for the closed universe. In 3c and 3d all the 
data are same with 3b except for 7 which has values respectively 0.92 and 3.51 in 
these two cases. 3e and 3f represents the graph of the tachyonic field potential B 
against cosmic time t at k=-l,i.e. open universe and 3g-3h represents the tachyonic 
field potential B at k=0, i.e. flat universe. 
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2.2.1 Normal Tachyonic Scalar field (e = +1) : 

Here for real if) both numerator and denominator in equation (|26[) must have the same 
sign. The denominator cannot be negative because then the expression within the second 
square root in equation (|27|l becomes negative and hence B(ip) becomes imaginary. So, the 
numerator has to be positive. Hence, when the perfect fluid matter component is not of 
exotic nature (i.e., does not violate the weak energy condition) then emergent scenario is 
not possible in this case for flat Bianchi I (k = 0)model and open Bianchi III (k = —1) 
model whereas for closed Kantowaski Sachs model(fc = +1) emergent scenario is possible 
provided inequality (|18p holds. When the perfect fluid matter component is of phantom 
nature then emergent scenario is possible for all three models provided restriction (I19|l 



2.2.2 Phantom Tachyonic Scalar Field (e = -1) : 

In this case, for real if) the numerator and the denominator of equation (|26|) should have 
opposite signs. But the denominator can not be negative as stated earlier. So, the only 



So, when the perfect fluid matter component is not of exotic nature then emergent scenario 
is always possible when k — 0, — 1. While for k = +1, emergent scenario is possible when 
restriction (|20[) holds. When the perfect fluid matter component is of phantom nature, 
then emergent scenario is possible for fiat, closed and open models provided restriction 
(2TJ holds. 

3 Lemaitre-Tolman-Bondi model 

By introducing the mass function F(r, t) (Joshi 1979; Banerjee 2003; Debnath 2003,2004) 
(related to the mass contained within the co-moving radius r)as 



holds. 



Also if) can be obtained as 




2 



dt 




F(r,t) = R(R 2 - f(r)) 



(28) 



the Einstein's equations in LTB model are 



8nGp 



F'(r,t) 
R 2 R' 



(29) 



and 



8nGp 



F(r,t) 
R 2 R 



(30) 



and the evolution equation for R is 



2.R.R + R 2 + 8irGpR 2 = f(r) 



(31) 
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The conservation equation is 

p + 3H(p + p) = (32) 

where 

H = l + (33) 



3 \R' R / 

is the Hubble parameter. The equation of state is given by 

p = ep, where e is a constant. (34) 
Using equations ([33|) and (J34J, equation ([32|) can be integrated to give 



P = Po 



(R 2 R')~ (1+C) (35) 



with po, a constant of integration. 

Using (|34[1 and (|35[1 in equation (I31|) . the differential equation for i? is obtained as 

2flj? + fi 2 + =/(r) (36) 

Now, taking 87rG = 1 and writing 

= h(r) fl (t) (37) 

we can write (|36p as 

m) (5W) 3(1+e) i(t)+(» 2 ( 9 (t)) 3(1+0 + — - - ep l - — - 4? (fW) 3(1+e) 

dt U * ' (h(r) f ((h(r)) 2 &h(r)y 1+c) h2 

(38) 

Choosing e = — 1 (i.e., matter as a cosmological constant) the above differential equation 
can be written in separable form as 

i.e., 

„ / x rf 2 / s / d " 2 
2g(t)—g(t) + 

and 



ZffW^flW+^sW) =A (39) 



(h(r)f 
Solving (|39|) we have 



' iPo + /(r)} = A (40) 



(VAffW-l + v/¥0) 5 exp I (Aff(t) - 5) ^ffM = (t — to)} (41) 



and 



Kr) = (42) 
where to appears as the constant of integration. From the solution (|41|l , we note that in 
the asymptotic region(i.e.t — > —00), R(t,r) has a finite non-zero value. (g(t) — > i). Also, 
the asymptotic behaviour of g(t) is shown in figure 4. 

Thus, emergent scenario is possible in LTB model with matter only in the form of a 
cosmological constant. 
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Fig 4 

4 shows the variation of the function g(t) w.r.t. time t. 
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4 Discussion 



In this work, first we have examined the possibility of emergent scenario for general hon- 
ogeneous and anisotropic model of the universe filled with matter having non-interacting 
two components -one in the form of perfect fluid with linear equation of state p = Lop and 
a scalar field (real, phantom or tachyonic) with potential as the other component. For 
real scalar field and tachyonic scalar field , when the perfect fluid matter component is not 
of exotic nature (i.e. does not violate the weak energy condition) then emergent scenario 
is not possible for flat Bianchi I model (k = 0) and open Bianchi III model (k = — 1) 
but for closed Kantowaski-Sachs model (k = +1), emergent scenario is possible when in- 
equality (|f 8|) holds. When the perfect fluid matter component is of phantom nature,then 
emergent scenario is possible for all models in real scalar field and tachyonic scalar field 
provided restriction (|19|) holds. For phantom scalar field and phantom tachyonic scalar 
field, when the perfect fluid matter component is not of exotic nature, then emergent 
scenario is always possible for flat and open models while for closed model it is possible 
when inequality (|20|) holds. But when the perfect fluid matter component is of phantom 
nature, then emergent scenario is possible for all three models provided restriction (|21|) 
holds. As these restrictions for emergent scenario depend on the anisotropy scalar a so 
the results are distinct from the isotropic model. 

In the second case we have studied the emergent scenario for LTB model. Here it is 
possible to have a solution assuming the scale factor to be inseparable (product form) and 
the matter is purely a cosmological constant. The asymptotic analysis shows that in the 
infinite past there is no singularity and the scale factor grows with time. In future it will 
be interesting to study the emergent scenario with general form of matter, particularly, 
with exotic matter. 
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